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ABSTRACT 

The existence of various anomalous stars, such as the first stars in the uni- 
verse or stars produced by stellar mergers, has been recently proposed. Some 
of these stars will result in black hole formation. In this study, we investigate 
iron core collapse and black hole formation systematically for the iron-core mass 
range of 3 - SOMq, which has not been studied well so far. Models used here 
are mostly isentropic iron cores that may be produced in merged stars in the 
present universe but we also employ a model that is meant for a Population III 
star and is obtained by evolutionary calculation. We solve numerically the gen- 
eral relativistic hydrodynamics and neutrino transfer equations simultaneously, 
treating neutrino reactions in detail under spherical symmetry. As a result, we 
find that massive iron cores with ~ lOM© unexpectedly produce a bounce owing 
to the thermal pressure of nucleons before black hole formation. The features of 
neutrino signals emitted from such massive iron cores differ in time evolution and 
spectrum from those of ordinary supernovae. Firstly, the neutronization burst is 
less remarkable or disappears completely for more massive models because the 
density is lower at the bounce. Secondly, the spectra of neutrinos, except the 
electron type, are softer owing to the electron-positron pair creation before the 
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bounce. We also study the effects of the initial density profile, finding that the 
larger the initial density gradient is, the more steeply the neutronization burst 
declines. Further more, we suggest a way to probe into the black hole progeni- 
tors from the neutrino emission and estimate the event number for the currently 
operating neutrino detectors. 

Subject headings: black hole physics — relativity — hydrodynamics — neutrinos 
— radiative transfer — methods: numerical 



1. Introduction 

Various anomalous stars, such as the first stars in the universe (so-called Population III 
stars) or stars produced by stellar mergers in stellar clusters, are being studied recently. As 
for the Population III stars, it is suggested theoretically that they are much more massive 
(M > IOOM0) than stars of later generations (e.g., Nakamura & Umemura 2001). On 
the other hand, A^-body simulations show that the runaway mergers of massive stars occur 
and that very massive (M > IOOMq) stars are formed in a young compact stellar cluster 
(e.g., Portegies Zwart 1999). Especially notable is a newly suggested formation scenario for 
supcrmassive black holes which requires the formation of intermediate-mass black holes by 
the collapse of merged stars in very compact stellar clusters (e.g., Ebisuzaki ct al. 2001). If 
these anomalous stars collapse to black holes without supernova explosions, it is supposedly 
difficult to be hard to probe into their progenitors. One possible means of such a probe is, we 
think, to examine the neutrinos emitted during the black hole formation. For this purpose, 
systematic studies on black hole formation, including effects of the neutrinos, are needed. 

So far, various numerical simulations of supernova explosions have been done by many 
authors. Similarly, numerical studies on black hole formation have also been recently pro- 
duced (e.g.. Fryer 1999; Linke et al. 2001; Fryer et al. 2001, hereafter FWHOl; Sekiguchi 
& Shibata 2005, hereafter SS05; Nakazato et al. 2006, hereafter NSY06; Sumiyoshi et 
al. 2006, hereafter SYSC06). Fryer (1999) classified core collapse into three types, a) 
Stars with M < 25Mq make explosions and produce neutron stars, b) Stars ranging 
25Mq ^ M < 4OM0 also result in explosions but produce black holes via fallback, c) 
For M > 40 M0, the shock produced at the bounce can neither propagate out of the core 
nor make explosions. In any case, the core bounces once. SYSC06 computed fully gen- 
eral relativistic hydrodynamics under spherical symmetry, taking into account the reactions 
and transports of neutrinos in detail and confirmed class c) for the collapse of a progenitor 
with 4OM0. On the other hand, much more massive stars result in black hole formation 
without bounce. SS05 studied the criterion for the collapse without bounce. Their com- 
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putations are fully general relativistic, and they investigated the dynamics systematically, 
varying the initial mass and rotation. They concluded that non-rotating iron cores with a 
mass of Miron ^ 2.2M0 collapse to black holes without bounce. However, they employed the 
phenomenological equation of state and did not consider the effects of neutrinos. 

There are studies also on a collapse of very massive stars in a context of the evolution 
of Population III stars. As mentioned already. Population III stars may be very massive, 
M > IOOMq. It is supposed that the pair creation of electrons and positrons makes a star 
unstable during the helium burning phase if they do not lose much of their mass during the 
quasi-static evolutions because of zero metaUicity. Stars with < 260Mq reverse the collapse 
by rapid nuclear burnings and explode to pieces, which are called pair-instability supernovae, 
while more massive stars cannot halt the collapse and form black holes (e.g., Heger et al. 
2003). Note that, however, these numbers are still uncertain at present (e.g., Ohkubo et al. 
2006). Assuming that Population III stars with M > 3OOM0 are formed and evolve without 
mass loss, FWHOl and NSY06 showed that they collapse without bounce for spherically 
symmetric models under fully general relativistic computations while NSY06 treated the 
neutrino transport more in detail than FWHOl. FWHOl also computed the collapse of a 
rotating star with 3OOM0 under Newtonian gravity and showed that it has a weak bounce 
and then recoUapses to a black hole immediately. As for the collapse of supermassive stars 
with M > 5 X IO^Mq, Linke et al. (2001) found that they form black holes without bounce 
before becoming opaque to neutrinos. 

The black hole formation of stars in the mass range between ~ lOOM© and ~ 26OM0, 
which corresponds to the iron-core mass range between ~ 3M0 and ~ 3OM0, has not 
been studied well so far. This is because they are supposed to explode as pair-instability 
supernovae during the quasi-static evolutions if they are single stars. Recently, on the other 
hand, stars produced by stellar mergers in a young compact stellar cluster were studied 
in detail and their evolutionary paths are beginning to be revealed (Suzuki et al. 2007). 
While they do not calculate the evolutions of these stars up to the black hole formation, we 
speculate, as in § 13.5. H that they may avoid the explosions as pair-instability supernovae 
and form a massive iron core of the above-mentioned range. Therefore, we investigate, in 
this study, the iron core collapses systematically for the iron core masses of 3M0 and 3OM0, 
although there is no evidence to show their existence so far. 

To be more specific, we assume that the mass of an iron core is mainly determined 
by the entropy per baryon, and our investigation is done systematically for entropy. We 
solve the general relativistic hydrodynamics under spherical symmetry. We also solve the 
neutrino transfer equations simultaneously, treating neutrino reactions in detail. In addition 
to the isentropic iron core models, we employ the realistic stellar model of lOOM© and zero 
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metallicity, supposedly a Population III star, by Nomoto et al. (2005). We address the issues 
concerning the black hole formation of the merged stars in connection with our study and 
estimate the neutrino event number for the currently operating detectors. We also suggest 
a way to probe into the progenitors from the detection. We hope that this study will be 
not only a reference for future multi-dimensional computations but also provide a basis for 
neutrino astrophysics in the black hole formation. 

2. Initial Models and Numerical Methods 

At first, we construct the iron core models, which will later be used as initial models 
for the dynamical simulation of the collapse. The progenitor with 4OM0 in SYSC06 has 
an entropy per baryon, s ~ l.dks and an iron core mass, Miron = 1.98M0, whereas the 
massive Population III star models in NSY06 have s > Wks and Miron > 50Mq. In this 
study, we intend to bridge the gap of the black hole progenitors and discuss the neutrino 
emission systematically for this range. Unfortunately, realistic models of the progenitors for 
this range are rare, and "systematic" models for them are absent up to the present, since 
their astrophysical counterparts are not well known, as mentioned already. Therefore, we 
construct the initial models by ourselves. 

We assume that the iron cores in equilibrium configurations collapse by photodisinte- 
gration, as is the case for the onset of ordinary core collapse supernovae. We obtain the 
initial models, solving the Oppenheimer-Volkoff equation with the equation of state by Shen 
et al. (1998a, 1998b) assuming isentropy and the electron fraction = 0.5 throughout the 
core. We define the mass of the iron core, Miron, as the mass coordinate where the temper- 
ature is 5 X 10^ K, whereas we set the outer boundary at a much larger radius so as not to 
affect the dynamics. For the systematic analysis, we set the initial central temperature as 
^Initial = 7.75 X 10^ K, which is slightly higher than the critical temperature for the photo- 
disintegration (Figured]), and generate models la-6a with the values of entropy per baryon, 
s = 3kB-13kB, which have not been studied well so far, as mentioned above. In order to 
investigate the ambiguity in the onset of collapse, we also adopt a model (model 2b) with the 
same initial entropy per baryon as model 2a (s = Aks) but having half the central density. 
The key parameters of these models are summarized in Table [H In addition, we also employ 
the realistic stellar model of IOOMq with a vanishing metallicity by Nomoto et al. (2005) in 
order to validate the isentropic models. This model is supposedly a Population III star and 
resides in the range s = Sks-lSkB. 

As a next step, we compute the dynamics of spherically symmetric gravitational collapse 
with the neutrino transport. As for our numerical methods, we follow NSY06 and use the 



- 5 - 



general relativistic implicit Lagrangian hydrodynamics code, which solves simultaneously 
the neutrino Boltzmann equations (Yamada 1997 ; Yamada et al. 1999 ; Sumiyoshi et al. 
2005). We consider four species of neutrino, i/g, ^'e, and u^, assuming that Uj. and 
are the same as and P^, respectively, and take into account 9 neutrino reactions listed in 
NSY06. We use 127 radial mesh points, while 12 and 4 mesh points are used for energy- 
and angular distribution of neutrino, respectively. In order to assess the convergence of our 
results, we compute models with higher resolutions. They have the same initial conditions 
as model 2a. For model 2m, the number of radial mesh points is increased to 255. Model 2e 
use 18 mesh points for the energy spectrum while model 2g has 6 mesh points for the angular 
distribution. 

It is noted that our method allows us to follow the dynamics with no difficulty up 
to the apparent horizon formation. The existence of the apparent horizon is the sufficient 
condition for the formation of a black hole (or, equivalently, of an event horizon). For the 
Misner-Sharp metric (Misner & Sharp 1964) adopted in our computations, the radius of the 
apparent horizon is written as 



where c and G are the velocity of light and the gravitational constant, respectively (van 
Riper 1979). r is the circumference radius and fh is the gravitational mass inside r. Since 
our models are spherically symmetric, there is no difficulty in finding the horizon. 



3. Results and Discussions 

In this section, we show the results of our computations and discuss them. We study 
the dynamics of the collapse in § 13.11 and investigate the features of the neutrinos emitted 
during the collapse in § 13.21 In § 13. 3[ we investigate the role of the initial velocity or the 
deviation from equilibrium. We also make a comparison with the realistic progenitor models 
in § 13. 4[ Finally, we mention about the astronomical counterparts of our models and the 
possibility of the probe into the progenitors of the events in § 13.51 

To overview the characteristics of the models which we surveyed, we show in Figuredl the 
evolution of the central density and temperature of our results together with those of other 
simulations of black hole formation. The trajectories of the current models shown by solid 
lines are between those of previous models reflecting the different values of entropy. From 
this flgure, we can recognize that our investigation bridges the gap between two previous 
studies, SYSC06 and NSY06. 
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3.1. Dynamical Features 

It is known tliat ordinary supernovae with s ~ I/cb bounce because their central density 
exceeds the nuclear density (~ 2.5 x lO^'^g cm~^) and pressure drastically increases. From our 
computations, we find that models with 3/cb < < T.Sfce {M < IO.GMq) have a bounce and 
that they recoUapse to black holes. On the other hand, models with s > T.bk^ (M > 10.6Mq) 
collapse to black holes directly without bounce. We show the evolution of core collapse in 
Figure H] for two representative cases. 

In the case of 3k^ < s < 7.5A;b, it is noted that the bounce mechanism of the core with 
s > S/cb is not the same as that of ordinary supernovae. The high entropy cores bounce 
because of the thermal pressure of nucleons at sub-nuclear density. We can see this fact 
from the evolutions of central density and temperature in the phase diagram of the nuclear 
matter at Ye = 0.4 and 0.2 (Figure [3]). We note that for all models at the center. Ye ~ 0.4 
and Ye ~ 0.2 when T ~ 1 MeV and T ~ 10 MeV, respectively. These figures show that 
the models with higher entropies go from the non-uniform mixed phase of nuclei and free 
nucleons to the classical ideal gas phase of thermal nucleons and a particles, whereas that of 
an ordinary supernova goes into the uniform nuclear matter phase. In the ideal gas phase, 
the number of non-relativistic nucleons and a particles is comparable to that of relativistic 
electrons. Since the adiabatic index of non-relativistic gas is 7 = | and that of relativistic 
gas is 7 = |, the collapse is halted and bounce occurs. 

Because this bounce is weak and the shock is stalled, the inner core (or the protoneutron 
star) grows beyond the maximum mass of the neutron star and recollapses to a black hole 
soon (left panel of Figure [2]). In Figure HI we show the maximum mass of the neutron star 
assuming isentropy and the constant electron fraction (Ye = 0.1) under the equation of state 
by Shen et al. (1998a, 1998b). It is noted that the maximum mass is larger than 3Mq for 
the neutron star with high entropies, s > 4A;b- Since the maximum mass of the neutron star 
depends on the equation of state, it should be remind that the time interval from the bounce 
to the recollapse also depends on it (SYSC06). We will refer to this point again later. 

In Table [H we show the inner core mass, central density, temperature and adiabatic 
index at the bounce together with the interval time from the bounce to the apparent horizon 
formation. We can recognize that the density and the adiabatic index at the bounce get lower 
for the models with higher initial entropies. These features indicate that the bounce is not 
due to the nuclear force but to the thermal pressure of non-relativistic gas for high entropy 
cores. Moreover, the interval time from the bounce to the apparent horizon formation is 
shorter for the higher entropy cores. This is because the initial mass of the iron core (Miron) 
is larger than the maximum mass of the neutron star (Mmax) for the models with high 
entropies and they can collapse to black holes quickly. 



- 7- 



We also show the results for the models with higher resolutions in Table [H The central 
density and the adiabatic index at bounce, which are key parameters in our analysis, are 
not very different for models 2a, 2g and 2m. The central density at bounce of model 2e, 
which has 1.5 times finer energy mesh, is different by 14% from that of model 2a. This is 
because neutrinos affect the entropy variations before the neutrino trapping. In fact, the 
central entropy at bounce of model 2a is 3.50A;b while that of model 2e is 3.62A;b- However, 
qualitative features of their bounces are not changed. On the other hand, the interval times 
from the bounce to the apparent horizon formation are different by < 15% for models 2a, 
2e, 2g and 2m. This is because the start point of the recoUpase is roughly determined by 
the maximum mass of the neutron star as mentioned already. Since the mass accretion 
rate is of the order of lOM© s~^ during this phase in our models, the difference of O.IM© in 
the maximum mass means the difference of 10ms in the interval time, which is close to the 
discrepancies found here. Thus, the precise determination of the interval time is difficult in 
general. However, its dependence on the initial entropy is well established. 

We compare our results with other studies. In SS05, the non-rotating models with 
an iron core mass of > 2.28Mq end up with black holes without bounce. In our models, 
on the other hand, it is shown that the iron core with < 10. GM© (or the initial entropy 
s < 7.5/cb) has bounce before black hole formation. This discrepancy comes from the fact 
that their equation of state is parametric and does not take into account properly the effects 
of thermal nucleons in the collapsing phase. On the other hand, the rotating Population HI 
star with ~ SOOMq has a weak bounce at ~ lO^^gcm"^ in FWHOl, and these authors 
adopt a realistic equation of state (Herant et al. 1994). Since the rotation tends to produce a 
bounce, we can predict that the bounce is inevitable for an iron core with the mass < lOM© 
irrespective of rotations, and the effects of thermal nucleons are crucial. 

In the high entropy case s > 7.5/cb, more massive cores do not have a bounce but form 
an accretion shock before the apparent horizon formation. This is because, the outer region 
keeps collapsing supersonically while the central region becomes gravitationally stable by 
the thermal pressure of non-relativistic gas. We can see this feature in the right panel of 
Figure [21 As the initial mass gets larger, the transition occurs smoothly from the collapse 
with bounce to the one without bounce. Incidentally, the features of direct collapse are 
almost the same as those for the Population HI models in NSY06, where a detailed analysis 
can be found. 
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3.2. Neutrino Signals 

In this section, we discuss neutrino emission during core collapse. As mentioned already, 
we compute the collapse until the formation of the apparent horizon. However, the location 
of the event horizon is not known for our models although it is proved mathematically that 
the event horizon is always located outside the apparent horizon. Moreover, the numerical 
difficulty prevents us from computing the dynamics until the apparent horizon swallows 
the shock surface entirely. Because of these facts, the total energy and number of emitted 
neutrinos have some ambiguities. In this study, we estimate the upper and lower limits for 
the total energy and number of emitted neutrinos, following NSY06. The upper limit is 
obtained with an assumption that all neutrinos in the region between the shock surface and 
the neutrino sphere flow out without being absorbed or scattered. For the lower limit, on 
the other hand, we assume that all neutrinos in this region are trapped and do not come 
out. Fortunately, for the models with bounce, these ambiguities are minor, compared with 
the direct collapse models in NSY06, since the duration from the shock formation to the 
apparent horizon formation is longer and almost all neutrinos are emitted during this phase. 

The calculated results of the neutrino emission are summarized in Table [2l It is noted 
that we assume that Vr {t^t) is the same as i^fi), and that the luminosities of and 
are almost identical because they have the same reactions and because the difference of 
coupling constants is minor. In the following, ignoring this tiny difference, we denote these 
four species as collectively. In Table [21 we can recognize that the total energy does not 
change monotonically with the initial entropy of the core. This is because the duration of 
the neutrino emission is longer for the lower entropy models, while the duration is shorter 
and the neutrino luminosity is larger for the higher entropy models. 

In Figure [5l we show the time evolutions of neutrino luminosity for several models 
under the assumption that the neutrinos outside the neutrino sphere flow freely after the 
apparent horizon formation. As already mentioned, the time interval from the bounce to 
the apparent horizon formation depends on the equation of state. In SYSC06, it is shown 
that the features of the neutrino emission, such as a neutronization burst, are not sensitive 
to the equation of state very much for the early phase. From Figure O we can see that the 
sign of neutronization burst becomes less remarkable and disappears for the higher entropy 
models. 

In order to analyze these features, we discuss the neutrino emission from model la, as a 
reference model. In the upper left panel of Figure [6], we show snapshots of the luminosity of 
an electron-type neutrino as a function of the baryon mass coordinate. We can recognize that 
neutrinos are emitted on the shock surface mainly. The luminosity on the shock surface has 
a peak (e.g., at 1.25Mq in the upper left panel of Figure E]), which is similar to the situation 



- 9 - 



for ordinary supernovae (e.g., Thompson et al. 2003). In the following, we estimate the 
value of the luminosity semi-analytically and compare it with the results of our numerical 
simulations. 

At first, the number density of neutrinos on the shock surface can be evaluated roughly 
by the equilibrium value, 

neq.(e)(ie oc r de, (2) 

^^P ( w) + 1 

where T and [ly are the temperature and the chemical potential of the electron-type neutrino 
in /9-equilibrium at the shock surface, respectively, and is the Boltzmann constant. Here 
the value iiy is defined as /i^, = /ig — (/in — A*p), where /ig, [i-n and [ip are the chemical potentials 
of electron, neutron and proton, respectively, and they are given in the equation of state by 
Shen et al. (1998a, 1998b). The number flux is estimated as cneq. (cos6'), where {cos9) is 
a mean value of the angular cosine over the neutrino angular distribution and c is the light 
velocity. In Figure [TJ we compare the results of our numerical computation with the number 
flux estimated above. We can see that the equilibrium is not achieved completely, but the 
fraction is rather constant, ~ 0.6. Therefore, the luminosity is well estimated by 

L{e)de = C Y^J-^^)/ de, (3) 
h^c- (exp (^) + l) 

where h is the Planck constant, r is the radius of the shock surface and C ~ 0.6. 

From equation ([3]), we can see that the luminosity is determined by r, fi,^, T and (cos 6*). 
According to our numerical computation, T ~ 1.5 MeV and {cos 6) ~ 0.5 do not change very 
much on the time scale of the neutronization burst. Thus, the luminosity is dictated mainly 
by r and fi^. Snapshots of the proflles of /it, are shown in the lower left panel of Figure El 
and we can see that /ii/ has a peak on the shock surface for the following reason. When 
matter accretes onto the shock, the baryon mass density and the electron number density 
rise, leading to the increase of fie and, as a result, fi^- Immediately thereafter, neutronization 
occurs and the value of (/z„— /ip) rises, which reduces fi^- We can recognize from Figure[6]that 
the peaks of fi^, and the luminosity are correlated. As for the time evolution, the luminosity 
on the shock surface is lower at the early phase because the shock radius is small. On the 
other hand, it is also lower at the late phase because fi,^ is lower. This is the reason why the 
luminosity on the shock surface has a peak. 

We now investigate model 4a, whose initial entropy is s = 7.5A;b- In the lower right 
panel of Figure [6], snapshots of the proflles of fi^ for model 4a are shown. We can see that 
the value of //jy at the shock surface is lower than that of model la at the early phase. This 
is because the baryon mass density on the shock surface of model 4a at the bounce is lower 
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than that of model la, as has been mentioned. Accordingly, the electron number density 
and fie are also lower for model 4a, and fiy does not rise so high. This is the main reason why 
the neutronization burst is not remarkable. It is noted, moreover, that the electron fraction. 
Ye, on the shock surface of model 4a is lower than that of model la. This is because nuclei 
do not exist and the nucleons are already neutronized on the shock surface. The absence 
of nuclei is consistent with the fact that the higher the initial entropy is, the earlier nuclei 
dissolve into nucleons, as explained by Figure [31 In addition, we can see that the luminosity 
of Ue rises monotonically. This is because the area of a shock surface increases whereas fi^ is 
almost unchanging. 

The results for the models with higher resolutions are shown in Figure [81 While the 
duration times of their neutrino emissions differ slightly among the models as mentioned 
already, the profiles of their neutronization bursts are not very different qualitatively. In 
fact, the luminosity declines a little after the peak and increases again for model 2a. This 
feature is well kept in other models with higher resolutions. 

We show the time-integrated neutrino spectra in Figure M We can see that the spectra 
become softer for higher entropy models, especially for Uf, and u^. In order to investigate 
this tendency, we show the time-integrated spectra of the neutrino emitted before and after 
the shock formation in Figure [TOl We can see that, for higher entropy models, Ue and 
are also emitted before the shock formation. They are created by the electron-positron pair 
annihilation, and their energy is relatively lower (< several MeV) because the temperature 
is low (T < 1 MeV). On the other hand, for lower entropy models, Ue and can not 
be produced by the electron-positron pair process because positrons are absent owing to 
Pauli blocking. As for the Ue and Ux emitted after the shock formation, they are mainly 
created by bremsstrahlung. In this phase, the temperature near the neutrino sphere rises to 
T ~ several MeV, which makes the neutrino energies relatively high: ~ 10 MeV. Since the 
low energy (< several MeV) neutrinos are not emitted to any great extent and the spectra 
become harder for lower entropy models, the emission of low energy z/g and is characteristic 
for the collapse of high entropy cores. 



3.3. Initial Velocity Dependence 

We compare the results of models 2a and 2b, which are different in the initial values 
of central density and temperature, but have the same initial values of entropy per baryon 
[s = 4/cb)- We can consider that models 2a and 2b are the same model but with different 
initial velocities, because the density profile of model 2b at the time when the central density 
reaches that of the initial model of 2a almost coincides with that of model 2a (Figure [TT]) . 
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In reality, the onset of a collapse is determined not only by the core structure but also by 
the whole stellar structure. Thus, studying the initial velocity dependence of the core is 
meaningful. 

As a result of this comparison, we find that the initial velocity does not affect crucially 
the ensuing dynamics and the features of emitted neutrinos such as total number spectra or 
the time evolutions of the luminosity. This is because the velocity of model 2b at the time 
in Figure [TT] is several times lower than the sound speed at each point. For instance, the 
fastest point of model 2b in Figure [TT] has the velocity ~ lO^cms"^ while the sound speed 
is ~ 7 X lO^cms"^, there. If the supersonic region, where the infalling velocity exceeds the 
sound speed, existed in the initial model, the initial velocity profile may be important for 
the dynamics. However, since the temperature of our initial models is slightly higher than 
the critical temperature for the photodisintegration instability, they are unlikely to have 
supersonic region. 



3.4. Collapse of Population III Star with lOOM© 

In this section, we consider yet another example of very massive stars, that is, a Pop- 
ulation III star with lOOM©. We use a model constructed by Nomoto et al. (2005) with 
evolutionary calculations, and we refer to it as model R. This model is very massive and its 
entropy at the center is higher than that of ordinary supernova progenitors when it starts 
to collapse because the star does not lose its mass at all in its evolution owing to its zero 
metallicity. It should be emphasized that the isentropic models are meant for the massive 
stars that may be produced in the present universe, for example, by stellar mergers in clus- 
ters whereas Model R corresponds to a first-generation star in the past universe. Here we 
are interested in the differences that these models may make. In Figure [121 we show the 
comparison of the initial state of model R and our isentropic models at the time when their 
central densities become the same as that of model R. We can recognize that model R has the 
entropy ~ S.Sfce in the central region, which is between those of model la and 2a, whereas 
the iron core of model R is smaller than that of our models. In fact, the iron core mass of 
model R is ~ 2.32M0, which is close to that of model la. We show some of the initial values 
at the center of model R in Table [H Incidentally, the initial velocity profile is taken into 
account for model R although it is much lower than the sound speed at each point. 

As a result of collapse, model R has a bounce and recollapses to a black hole. As shown 
in Table [H the values of the central density and the central adiabatic index of model R at the 
bounce are between those of model la and 2a. This suggests that these values are determined 
by the initial central entropy as mentioned in § 13. 1[ On the other hand, model R has a much 
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longer time interval from the bounce to the recollapse, compared with our models. This is 
because the inner core mass of model R at the bounce (Mbounce) is smaller and the lower 
density of the outer core (Figure [T2D gives lower accretion rates. For instance, at t = 0.06 s, 
model R has a mass accretion rate ~ s~^ at the shock surface whereas model la has 
~ IIM0 s~^. Thus it takes much time until the inner core mass exceeds the maximum mass 
of the neutron star. 

We show the total energy of neutrinos emitted during the collapse of model R in Tabled 
and the time evolution of the emitted neutrino luminosity in Figure [121 We can see that 
the total energy of emitted neutrinos is larger than that of other isentropic models, despite 
the fact that the neutrino luminosity of model R is lower than those of our models. This 
is because model R neutrino emission lasts much longer. Moreover, the mean energy of 
the emitted neutrinos is larger for model R. This is also due to the longer duration time. 
The neutrino spectrum gets harder in the late phase because the density of the accreting 
matter becomes lower and the temperature on the neutrino sphere gets higher. Thus, the 
longer the duration time of neutrino emission is, the larger the mean energy of the emitted 
neutrinos becomes. It is noted that the duration time is sensitive to the equation of state, 
which is already mentioned, and hence the total and mean energy of emitted neutrinos is 
also sensitive to the equation of state. 

In the following, we discuss the features of the emitted neutrinos from model R for the 
early phase, which is not sensitive to the equation of state as already mentioned. Comparing 
Figures O and [121 we can see that for model R, the peak luminosity of the electron-type 
neutrino by the neutronization burst is lower than those of our models. The reason why it is 
lower than that of model la (s = 3A;b) is because the chemical potential of an electron-type 
neutrino for model R is lower than that for model la, while the shock radii in both models 
are not so different from each other (right panels of Figure [T^ . It is consistent with the fact 
that the density at the bounce of model R is lower than that of model la (Table [T]). On 
the other hand, the shock radii of models with s > Ak-B (models 2a-4a) are larger than that 
of model R. This is the reason why the luminosity of the neutronization burst for model R 
is lower than those of models 2a-4a. Furthermore, since the outer core density of model R 
is much lower than those of isentropic models (Figure [T^ . /j.^, drops quickly and the shock 
radius does not get much larger after the neutronization. From equation ([3D, These features 
lead to the fact that the luminosity of the electron-type neutrino after the neutronization 
burst drops more steeply for model R than for our models. It follows, then, that the decline 
of the neutronization burst depends not only on the initial entropy but also on the initial 
density profile. In particular, the larger the initial density gradient is, the more steeply the 
neutronization burst declines. 
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To sum up, the key parameters listed in Table [T] at bounce (e.g. central density, tem- 
perature etc.) do not differ very much between the isentropic models and model R. This 
is not true for the time profile of the neutronization burst because they depend not only 
on the central density at bounce but also on the initial density profile. However, since, in 
general, more massive iron cores have larger entropies, the following trend is generally true: 
The neutronization burst will become less remarkable as the progenitor gets more massive. 



3.5. Astrophysical Implications 

3.5.1. Progenitor of IMBH 

For supermassive black holes (SMBH) located at the center of many galaxies including 
ours, a new formation scenario via intermediate-mass black holes (IMBH) has recently been 
suggested (e.g., Ebisuzaki et al. 2001, Portegies Zwart et al. 2006). According to this 
scenario, very dense stellar clusters are initially formed in the vicinity of the galactic center 
(< 10 pc), and the massive stars with ~ 2OM0 in them undergo runaway collisions to form 
IMBHs before they lose most of their mass by supernova explosions and/or pulsations. After 
that, these IMBHs merge together and finally form SMBH. This scenario is supported by 
the discovery of the ultra luminous X-ray compact sources in M82 galaxy, which indicate the 
existence of IMBHs. It is conceivable that similar events occur in the Milky Way Galaxy as 
mentioned later. This scenario assumes that the supermassive stars formed by the runaway 
collisions would collapse to IMBHs when they are ~ IOOOMq. 

Recently Suzuki et al. (2007) have studied the structures and evolutions of these merged 
stars in the hydrogen burning. According to them, the smaller star sits at the center of the 
larger star after the merger of two stars with different masses. It is also demonstrated that the 
merged stars become convectively unstable by the positive gradient of the mean molecular 
weight and that their evolutions thereafter approach those of the single homogeneous star 
with the same mass and abundance. The central entropies of the merged stars will then 
be larger than those of the inhomogeneous single stars with the same mass. This suggests 
the possibility to form the IMBH progenitors by the merger without experiencing the pair 
instability. Here we speculate the entropy of these stars using previous studies on the single 
Population HI stars. Since the iron core of the Population HI star with lOOM© has entropy 
of ~ 3.5A;b (Nomoto et al. 2005), it is expected that these IMBH progenitors have entropies 
>3A;b. 

Even if the pair instability occurs, the massive stars corresponding to our models may 
still be formed. In fact, the positive entropy gradient and/or rotation may suppress the 
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convection in the merged star and the entropy at the center may remain low after the merge. 
Then the merged star has a massive envelope with a smaller core than the single stars with 
the same total mass. If the pair instability occurs for these objects, the nuclear burning 
may not produce total disruptions but lead to the eventual collapse. Again inferring from 
single Population 111 stars, we speculate that the central entropies of the merged stars will be 
smaller than s ~ IQks, which corresponds to SOOM© in NSY06. It is incidentally mentioned 
that the relations between the total mass and the iron core mass of merged IMBH progenitors 
is highly uncertain at present. 

In the preceding sections, we have shown that the neutrino signals from the black hole 
formation arc sensitive to the inner region of the progenitor. In this section, assuming our 
models correspond to above-mentioned merged stars which collapse to IMBHs at the center 
of our Galaxy (~ 8.5 kpc from the sun), we estimate the neutrino event number for the 
currently operating detectors. 

As for the event rate of the IMBH formation, based on the above-mentioned scenario 
and the fact that the SMBH residing in the center of our Galaxy (SgrA*) is ~ 3.5 x IO^Mq 
and the age of our Galaxy is ~ 10 Gyr, a very rough estimation for the formation rate of 
IMBH with ~ IOOOM0 is < once per 1 Myr. It is, however, mentioned that this event rate 
may be underestimated because star formation may not be continuous but triggered by some 
environmental effects (e.g., the merger of galaxies). Recent observations by Paumard et al. 
(2006) have revealed the existence of about 80 young massive stars within a distance of a 
parscc from SgrA* and some of them arc identified as OB stars and their ages arc about 6 ± 2 
Myr. These facts indicate that stars are actively formed in this region at present. Moreover, 
the IMBH candidate with ~ ISOOM©, IRS 13, is found in the same region (Maillard et al. 
2004). Thus, SgrA* may be currently growing under this scenario. 

In the following estimations for the neutrino event number, we do not take into account 
the neutrino mixing, although it should be. Since the mixing occurs mainly in the resonance 
regions and they are located outside the iron core of the progenitor, the neutrino oscillation 
does not affect the dynamics of core. Unfortunately the structures of the envelopes of 
merged stars, which are crucial for the neutrino mixing, are quite uncertain. There remain 
uncertainties as well on the mixing parameters, such as the mixing angle of sin^ 2^13 or the 
mass hierarchy. Thus, the precise evaluation of the neutrino flux including the neutrino 
mixing is deferred to future study. 
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3.5.2. Detection of low energy z/g hy Super Kamiokande and KamLAND 

As already mentioned, a good deal of low energy z/g is emitted from the collapse of 
the high entropy cores, which softens the spectrum. We estimate the Ve event number for 
Super Kamiokande III and KamLAND, currently operating neutrino detectors, under the 
assumption that the black hole formations considered in former sections occur at the center 
of our Galaxy. For both detectors, the dominant reaction is the inverse beta decay, 

Ue+V — > e"^ + n, (4) 

which we take into account only. We adopt the cross section for this reaction from Vogel & 
Beacom (1999). For Super Kamiokande III, we assume that the fiducial volume is 22.5 kton 
and the trigger efficiency is 100% at 4.5 MeV and 0% at 2.9 MeV, which are the values at 
the end of Super Kamiokande I (Hosaka et al. 2006). For KamLAND, we assume 1 kton 
fiducial mass, which means that 8.48 x 10^^ free protons are contained (Eguchi et al. 2003). 
We also assume that the trigger efficiency is 100% for all Ve energy larger than the threshold 
energy of the reaction. 

The results are given in Table [31 The total event number does not change monotonically 
with the initial entropy of the core because the total number of neutrinos depends on both 
the core mass and the duration time of neutrino emission, as already mentioned. In order 
to investigate the hardness of Ve spectrum, we calculate the ratio of the event number by Pe 
with < 10 MeV to that for all events. The ambiguity about the distance of source is also 
canceled by this normalization. This ratio gets larger as the entropy of the core becomes 
higher. This suggests that we can probe the entropy of the black hole progenitor especially 
in higher regimes [s > 7.5A;b) because the event numbers of Ue with < 10 MeV are over 100 
by Super Kamiokande III. 



3. 5. 3. Detection of neutronization burst by SNO 

The SNO detector consists of 1 kton of pure heavy water (D2O) and can distinguish ly^ 
flux by the charged-current reaction of the deuterium disintegration. Since SNO can also 
detect the Ue flux, we can estimate the intensity of the neutronization burst by comparing 
the event from the charged-current reaction of u^, 

Ue + d — >p + p + e~, (5) 

and that of Ug, 

Ue + d — > n + n + e'^ , (6) 
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using the SNO detector. SNO can also detect the neutral-current reaction, 

1/ + d — > n + p + u, (7) 

for all species. It is noted that the neutral-current reaction contains (= ly^, u^, Vj- and Ur) 
and the neutrino sphere of differs more from that of than that of z/g in general. Thus, 
for the comparison with reaction ([5]), reaction ([6]) is more appropriate than reaction ([7]). On 
the other hand, we also use ([7]) for the comparison because the event number of ([7]) is larger 
than that of ([HD. In our calculation, we use the cross sections from Ying et al. (1989) and 
assume that the trigger efficiency of these reactions is 100%. In fact, it is ~ 92% these days, 
which is the neutron (in the right hand side of equation ([7])) capture efficiency on ^^Cl and 
deuterons (Oser 2005). 

In the following analysis, we regard the emission of neutrinos before t = 0.06 s as the 
neutronization burst, where the time t is measured from the bounce. The criterion t = 0.06 s 
is chosen empirically from our simulations as an expedient. The method for extracting the 
neutronization burst from detection should be reconsidered for more detailed studies. Here 
we calculate the event numbers for t < 0.06 s as well as those for the entire duration time of 
the neutrino emission, and the results are summarized in Table HI We can recognize that the 
ratios of the i^e event number (A/",^^ <o.o6s) to the total event number of the charged-current 
reactions {N^^^<oms + ^Pe,<o.o6s) and that for the neutral-current reaction (A0vc,<o.o6s) are 
larger for the models whose neutronization burst declines more steeply. Despite the fact that 
these neutronization burst numbers are of the order of 10, we can probe into the black hole 
progenitors in principle. 

It is finally noted that the estimations in the current study are based on the spherically 
symmetric models. If progenitors are rotating rapidly, the neutrino sphere will become non- 
spherical and the neutrino emissions will be affected in general. This will be the subject of 
future investigations. 

4. Conclusions 

In this paper, we have numerically studied gravitational collapse and black hole forma- 
tion of massive iron cores systematically, taking into account the reactions and transports 
of neutrinos in detail. Massive iron cores with ~ IOMq have a bounce owing to thermal 
nucleons, following which they collapse to black holes when the maximum mass is reached. 
As for the emitted neutrinos, the spectra of iJe and (= f/^, t'/^, t^t and Pr) become softer 
for more massive models, or higher entropy models, because a high entropy generates a large 
number of electron-positron pairs, which create Ue and u^- The neutronization burst from 
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more massive iron core becomes less remarkable or disappears completely. This is because 
the density at the bounce is lower and even the v^, number density in equilibrium becomes 
lower. 

We have found that if the initial velocity is lower than the sound speed, it does not 
affect the collapse very much. We have also compared the collapse of our isentropic models 
with that of the realistic model, which is obtained by the detailed modeling of the evolution 
of Population III stars and we have found that the steep decline of the neutronization burst 
depends not only on the initial entropy but also on the initial density profile. Moreover, 
assuming our models as the progenitors of IMBHs collapsing at the Galactic center, we have 
estimated the neutrino event numbers. As a result, for Super Kamiokande III, the ratio of 
the Pe event number for < 10 MeV to that for all events gets larger as the entropy of the core 
becomes higher, especially for s > 7.5A;b- We have suggested that we can use these features 
to probe into the progenitors. As for the lower entropy cores, despite the fact that the event 
number for the early phase of the emission is less than 100 by SNO, we have suggested that 
the steep decline of the neutronization burst can be distinguished in principle. 

Concerning the prediction of neutrino event number, there is a room for further improve- 
ment. Firstly, the effects of the neutrino oscillation should be taken into account. Secondly, 
multi-dimensional effects, such as rotation or magnetic field may be important, since they 
will affect the dynamics of collapse itself. This study will be hopefully prove a first step 
toward a neutrino astrophysics for black holes. 
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of the National Astronomical Observatory of Japan (VPP5000 System projects wknlOb, 
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project 108). This work was partially supported by Japan Society for Promotion of Science 
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Education, Science and Culture of Japan through 14740166, 15540243, 15740160, 17540267, 
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IQ4 IQ5 IQ6 IQ7 IQ8 IQ9 IQlO lO^^ 10'3 lO^^ lO^^ IQl^ 

density (g/cm^) 

Fig. 1. — The evolution of the central density and temperature for various models which 
result in black hole formation. The dot-dashed line is for a realistic progenitor with the 
initial mass AOMq in SYSC06 and the long-dashed line is for a Population III star with the 
initial mass IO5OOM0 (s = 74.75A;b) in NSY06. Three solid lines are for our models in this 
paper, and each line corresponds to model la (s = 3A;b), 3a (s = 5A;b) and 5a (s = IO/cb), 
from right to left. The shaded area represents a gravitationally unstable region by labeled 
physical processes. 



- 21 - 



Table 1: Key Parameters for all Models. 





'''initial 


-^^iron 


Pinitial 


^initial 


-^^bounce 


Pbouncc 


-^bounce 


Tbounce 


Precollapse 


model 


(fe) 


(Mo) 


(g cm~^) 


(K) 


(Me) 


(gcm^'^) 


(MeV) 




(msec) 


la 


3.0 


2.44 


2.71x10^ 


7.75x10^ 


0.75 


1.95x10^^ 


25.9 


2.38 


96.7 


2a 


4.0 


3.49 


1.40x10^ 


7.75x10^ 


1.10 


9.58x10^3 


26.9 


1.89 


62.0 


2b 


4.0 


2.93 


7.00x10^ 


6.86x10^ 


1.05 


9.90x10^3 


26.7 


1.91 


63.4 


3a 


5.0 


4.97 


8.82x10^ 


7.75 xlO^ 


1.5 


2.97x10^3 


19.4 


1.58 


52.6 


4a 


7.5 


10.6 


4.20x10^ 


7.75 xlO^ 


2.7 


3.00x10^2 


12.0 


1.54 


37.9 


5a 


10.0 


19.3 


2.67x10^ 


7.75 xlO^ 












6a 


13.0 


34.0 


1.84x10^ 


7.75x10^ 












R 


3.5 


2.32 


2.34x10^° 


1.61x10^° 


0.65 


1.37x10^-^ 


22.7 


2.19 


402 


2a 


4.0 


3.49 


1.40x10^ 


7.75x10^ 


1.10 


9.58x10^3 


26.9 


1.89 


62.0 


2m 


4.0 


3.49 


1.40x10^ 


7.75x10^ 


1.04 


9.64x10^3 


27.5 


1.88 


69.5 


2g 


4.0 


3.49 


1.40x10^ 


7.75x10^ 


1.05 


9.66x10^3 


27.2 


1.88 


65.7 


2e 


4.0 


3.49 


1.40x10^ 


7.75x10^ 


1.10 


8.25x10^3 


25.9 


1.80 


73.2 



Note. — Sinitiai IS the initial value of the entropy par baryon. Miron and Mbounce are the mass of initial 
iron core and inner core at the bounce {t = 0), respectively, pinitiai and Pbounce are the central density of the 
initial model and at the bounce, respectively. Ilnitiai and Tbounce are the central temperature of the initial 
model and at the bounce, respectively. Tbounce is the central adiabatic index at the bounce. trecoUapse is the 
interval time from the bounce to the apparent horizon formation. 

Table 2. Estimates of Average and Total Energies of Emitted Neutrinos. 



model 


{E.J 
(McV) 


(MeV) 


(MeV) 


TT'tot 

(lO^^crgs) 


TT'tot 

Ei?,,52 
(lO^^crgs) 


TT'tot 

(10-'''2crgs) 


ptot 
-^all,52 

(lO^^ergs) 


la 


11.01 - 11.01 


15.03 - 15.03 


19.60 - 19.60 


3.29 - 


3.29 


1.94 - 


1.94 


1.43 - 1.43 


10.96 - 


10.96 


2a 


10.32 - 10.32 


14.29 - 14.30 


19.30 - 19.30 


3.21 - 


3.21 


1.58 - 


1.58 


1.33 - 1.33 


10.12 - 


10.12 


2b 


10.17 - 10.17 


14.27 - 14.27 


19.35 - 19.35 


3.24 - 


3.24 


1.63 - 


1.63 


1.36 - 1.36 


10.29 - 


10.29 


3a 


9.29 - 9.29 


13.79 - 13.79 


19.31 - 19.31 


3.10 - 


3.10 


1.48 - 


1.48 


1.30 - 1.30 


9.78 - 


9.78 


4a 


7.30 - 7.30 


11.95 - 11.95 


19.55 - 19.55 


3.19 - 


3.19 


1.56 - 


1.56 


1.35 - 1.36 


10.18 - 


10.21 


5a 


6.24 - 6.25 


10.34 - 10.37 


18.37 - 18.70 


4.01 - 


4.02 


2.34 - 


2.35 


1.69 - 1.73 


13.11 - 


13.31 


6a 


5.24 - 5.25 


8.14 - 8.19 


14.33 - 14.33 


6.15 - 


6.17 


4.71 - 


4.75 


1.70 - 1.73 


17.66 - 


17.84 


R 


15.34 - 15.34 


18.90 - 18.90 


23.42 - 23.42 


9.42 - 


9.42 


7.89 - 


7.89 


4.40 - 4.40 


34.89 - 


34.90 



Note. — The mean energy of emitted Ui (with upper and lower limits) is denoted as (E.^) = E^J"^ /N^°^ , where 
and ArJ°* axe the total energy and number of neutrinos. is the total energy summed over all species. 
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time {s) time (s) 



Fig. 2. — Radial trajectories of mass elements. The left panel is for model 2a (s = Aks)', time 
is measured from the bounce. The right panel is for model 5a {s = lO/c^); time is measured 
from the point at which the apparent horizon is formed. 



Table 3: Event Numbers of from Super Kamiokande III and KamLAND. 



model 


J^i?e<W MoV,SK 




-/Viye<lOMcV,Kam 






^ Vf , Kam 


la 


3.3% 


6163 


3.3% 


174 


2a 


4.0% 


4778 


4.0% 


135 


2b 


4.0% 


4910 


4.0% 


139 


3a 


4.6% 


4319 


4.6% 


122 


4a 


7.3% 


4018 


7.3% 


114 


5a 


11.8% 


5326 


12.0% 


151 


6a 


20.1% 


9139 


20.5% 


259 



Note. — The subscript "< 10 MeV" means the event of with < 10 MeV, and the subscript "SK" and 
"Kam" mean the prediction for Super Kamiokande III and KamLAND, respectively. 
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log density (g/cm^) 



Fig. 3. — Phase diagram in p—T plane from Shen et al. (1998b) for fixed electron fraction, 
(thick lines). The nucleus exists in the region below these thick lines. The phase boundaries 
depend on Ye, whereas the same trajectories are plotted for the upper panel and the lower 
panel. The dashed line represents the evolution of the central density and temperature for 
the ordinary supernova progenitor with the initial mass 15Mq (Sumiyoshi et al. 2005), and 
the solid lines do the same for the progenitors studied. Each line corresponds to models la 
(s = Sk^), 3a (s = S^b) and 5a (s = IO/cb), from right to left. 
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2 4 6 8 

entropy (kg) 

Fig. 4. — The maximum mass of the neutron star assuming isentropy and a constant electron 
fraction which is isentropic (Ye = 0.1) under the equation of state by Shen et al. (1998a, 
1998b). The sohd and dashed hues represent the maximum mass in the sense of the baryon 
rest mass and the gravitational mass, respectively. 
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Fig. 5. — Luminosities of Ve (short-dashed hne), Ve (sohd hne) and Vx (long-dashed hne) as a 
function of t, where stands for /i- and r-neutrinos and their anti-particles. Squares show 
the time when the apparent horizon is formed. Upper left, upper right, lower left and lower 
right panels are for models la (s = 3A;b), 2a (s = 4A;b), 3a (s = 5A;b) and 4a (s = 7.5A;b), 
respectively. 
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Fig. 6. — Snapshots of the profiles for the luminosity and the chemical potential of an 
electron-type neutrino. The left panel corresponds to the model la (s = Sfce) and the right 
to the model 4a (s = 7.5A;b)- 



Table 4. Event Numbers by SNO. 



model 




.<0.06 s + Npe.<0 06 s 


-'VnC,<0.06 s 


.^NC,<0.06 a 






la 


30.7 6.4 


82.7% 


45.4 


67.5% 


84.2 


45.4 


201 


2a 


40.8 10.7 


79.2% 


70.4 


57.9% 


69.8 


30.1 


162 


2b 


42.6 12.4 


77.5% 


77.1 


55.3% 


73.6 


33.8 


179 


3a 


38.4 12.6 


75.3% 


78.5 


49.0% 


56.7 


24.7 


142 


4a 


31.8 15.1 


67.9% 


91.2 


34.9% 


37.4 


19.3 


124 


5a 










44.9 


32.7 


220 


6a 










59.3 


54.0 


221 


Note. 


— These values are the event 


number for the charj 


;ed-current reaction except AfNC,<0.06a 


and A^NC.all- 


The 



subscript "< 0.06 s" means the event at t < 0.06 s, where t is the time measured from the bounce, and the subscript "all" 
means the event for all duration times of neutrino emission. 
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Fig. 7. — Time evolutions of the number flux for the electron-type neutrino with the energy 
10 MeV < E < 20 MeV detected by the comoving observer. Solid lines and dashed lines 
represent the results of our computation and the values estimated by the number density 
in equihbrium and {cos 9) of our computation, respectively. The left panel corresponds to 
model la (s = 3/cb) at M = I.SMq and the right to model 4a (s = 7.5A;b) at M = 3.3Mq. 
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Fig. 8. — Luminosities of Ve as a function of t for models 2a (solid line), 2m (short-dashed 
line), 2e (long-dashed line) and 2g (dot-dashed line). The meaning of squares is the same as 
in Figure [H 
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Fig. 9. — Spectra of time-integrated emissions of v^, (short-dashed hne), v^, (sohd hne) and 
(long-dashed hne). Upper left, upper right, lower left and lower right panels are for models 2a 
(s = 4A;b), 4a (s = 7.5A;b), 5a (s = IO/cb) and 6a (s = 13/cb), respectively. 
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Fig. 10. — Spectra of time-integrated emissions of (short-dashed hne), Ue (sohd hne) and 
(long-dashed hne). The upper left and upper right panels give the time integrations of 
the emission before and after bounce, respectively, for model 2a {s = 4A;b). The lower left 
and lower right panels present the emission before and after shock formation, respectively, 
for model 6a (s = 13A;b)- 
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Fig. 11. — Comparisons of the density profiles (left) and the velocity profiles (right). Solid 
lines represent the initial profiles for model 2a; dashed lines represent the profiles for model 2b 
at the time when the central density becomes the same as the initial central density of 
model 2a. 
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Fig. 12. — Comparisons of the density profiles (left) and the entropy profiles (right). Solid 
lines represent the initial profiles for model R and dashed and dot-dashed lines represent the 
profiles for models la and 2a, respectively, at the time when the central density becomes the 
same as the initial central density of model R. 
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Fig. 13. — Results of the collapse for model R. In the left and right panel, the notations of 
lines are the same as in Figures O and El respectively, but the end points of the lines in the 
left panel represent the time when the apparent horizon is formed. 



